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SUMMARY 


Minimum weight structural design consists in determining the set of 
design parameters of the structure which makes the weight of the structure a 
minimum without violating the constraints imposed on the structure. 

In a structural optimization problem, where the constraints are the 
limiting stresses, failure is absolutely prevented, that is, with a probabil- 
ity of 1. Such absolute safety not only limits the minimization of weight but 
is idealistic and therefore impractical, for in actual systems limits must be 
imposed within which failure can occur. The limiting stress constraint, 
having a failure probability of 0, is hence replaced by a limiting probability 
of failure constraint. 

In this paper, the concept of probability of failure is explained for a 
general structural system under simultaneous and alternative load systems. It 
is then shown how the probability of failure constraint can be formulated in 
terms of the structural design parameters. This constraint is developed in 
detail for optimizing a simple truss whose design parameters are the section 
areas of the members. The optimizing weight function and the constraint in 
terms of the design parameters define the mathematical model of the structural 
optimization problem. The structural synthesis method is suggested for 
solving the mathematical programming problem, and it is shown how the proba- 
bility of failure constraint can be incorporated in the synthesis technique. 
The application of the constraint is illustrated by means of two simple 
examples with reference to the truss. 


INTRODUCTION 


Minimum weight structural design consists in determining the set of 
design parameters of the structure which makes the weight of the struct lore a 
minimum and at the same time ensures that the structure is safe under the 
governing load system. To determine whether a structure is safe or not, the 
computed design stresses are compared with the permissible stresses. The 
structure will fail if the stress in any component exceeds the allowable 
limit. 

In a structural optimization problem, where the constraints are the 
limiting stresses, the extent to which the weight can be diminished is 
governed, then, by the limiting stress constraint. Imposing this constraint 
prevents failure absolutely, that is, with a probability of 1. Now the design 




stresses in the components of a structure are a function of the loading, the 
structural design parameters, namely, the sizes of the elements and the 
material properties. The permissible stress or the failure stress for a 
particular material is a physical property of the material. The properties 
have to he estimated to represent all possible relevant conditions. Past 
experience shows that it is possible to represent a material property by a 
frequency distribution. Such a distribution, then, furnishes its most probable 
value, that is, its mean and its probable range of fluctuation represented by 
its standard deviation. Similarly, the applied loading is a random variable 
and may also be represented by a frequency distribution. This enables an 
assessment to be made of the most probable conditions, fluctuation ranges, and 
the probability that certain extremes will occur. 

It is seen, then, that to prevent failure absolutely, extreme fluctua- 
tions of load and permissible stresses would have to be taken into account, 
that is, the upper and lower extremes of the load and permissible stress. 

Such absolute safety not only limits the minimization of weight but is ideal- 
istic. On actual systems, limits must be imposed within which failure can 
occur. For the most conservative system, failure can be made highly improb- 
able by making its probability small but it is impractical to attempt to 
prevent it absolutely. The limiting stress constraint, having a failure 
probability of 0, is hence replaced by a limiting probability of failure 
constraint . 

In assigning probability of failure to a structure, a compromise has to 
be made between the level of economy to be maintained and the amount of risk 
a designer is prepared to undertake. Consider now, as examples, two struc- 
tural systems: (l) a system of limited purpose mission with a large number of 
units per mission, such as missiles, (2) a multipurpose single trial system, 
such as a Mars probe lander. In the case of missiles, say if 50 missiles were 
to be deposited over an area, it does not matter if 1 or 2 do not come off and 
fail to accomplish the mission. So it pays not to overdesign the missile and 
to allow a possibility of failure of 1 or 2 out of 50 as a result of stresses 
exceeding their limiting values. This means that in the design of the missile, 
the probability of failure should be less than or equal to 0.02 or 0.04. Such 
a constraint, as opposed to a failure probability of 0, results in consider- 
able saving in the weight of the missile. In the case of a single trial 
system such as a Mars probe lander, the probability of failure would have to 
be extremely small, say 1 in 10 7 . 

A reliability approach to structural design has been recognized increas- 
ingly during the past decade (see Bibliography); the researchers have stressed 
the problems of introducing realistic safety factors and of proportioning a 
structure to meet a prescribed reliability. This literature indicates the 
need for a method of systematic determination of the minimum weight design ‘of 
a structure within limits of a prescribed probability of failure. This paper, 
then, endeavors to present a compact mathematical model of the structural 
optimization problem wherein the merit function to be optimized is the struc- 
tural weight, and the probability of failure of the structure is expressed in 
terms of the design parameters and employed as a constraint in the 
optimization problem. 
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In this paper, me concept of probability of failure is explained for a 
general structural system under simultaneous and alternative load systems. 

It is then shown how the probability of failure constraint can be formulated 
in terms of the structural design parameters. This constraint is developed in 
detail for the optimization of a simple truss whose design parameters are the 
section areas of the members. The mathematical model of the structural opti- 
mization problem Is defined by the optimizing weight function and the con- 
straint in terms of the design parameters. The structural synthesis method is 
suggested for solving the mathematical programming problem, and it is shown 
how the probability of failure constraint can be incorporated in the synthesis 
technique. The application of the constraint is illustrated by means of two 
simple examples with reference to the truss. 


NOTATION 


A± 

E 


Qj(J=1,2) 


g u 


l 



ith design parameter (area of the ith member of the truss, with 
i =1, 2 for bars OA and OB, respectively) 

modulus of elasticity 

failure stress of the ith element, taken to be a normally 
distributed random variable 

probability of failure of the ith element under qth loading 
probability of failure of the structure under loading Q. 

probability of failure of the structure 

jth loading, taken to be a normally distributed random variable 

mean of random variable Q. 

J 

force in ith member under jth loading 

structural weight 

applied stress in ith element 

random variable [F - ( R l ] 

length of either member of the truss 

density of material for the truss 

standard deviation of random variable Q . 

J 
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ANALYSIS 


Factors Influencing the Probability of 
Failure of the Structure 

Dead load * - Factors causing variation in dead load are the dimensions of 
the structure and the specific -weight of the materials. 

Live load . - To have a complete representation of live load, it is neces- 
sary to know the most probable conditions of loading -with a high frequency of 
occurrence and, at the same time, to allow for certain adverse conditions that 
could occur. For example, if the loading represents wind load, one must know 
what percentage of time the wind velocity is in a certain range and also the 
magnitude of, say, hundred -year gales. 

The variation in loading, its most probable value, and estimates of cer- 
tain extreme cases can be represented by a frequency distribution with given 
mean and standard deviation. 

Structural stiff ness.- The structural stiffness depends upon the struc- 
tural configuration and the elasticity of the members. The modulus of elas- 
ticity of the structure is a constant, statistically only. It can be 
represented by a frequency distribution with given mean and standard devia- 
tion. The sizes and the dimensions of the structures can have definite 
variations about their mean values if their tolerances are specified. The 
design stress in a member is then a function of the loading and the 
structural stiffness. 

Failure stres s . - The failure stress is a physical property of the 
material. Its variability depends on the control exercised during 
manufacturing process. It, too, is represented by a frequency distribution. 

Development of Probability of Failure Constraint 

Consider a structure made up of r elements or components. Let Q1 
and Q2 represent two independent loading systems; Q1 and Q2 are random 
variables, normally distributed with means Q1 and Q2 and standard deviations 
cjq^ and 0 q 2> respectively. It could happen that both Q1 and Q2 would act 
simultaneously as dead and live loads on the structure. On the other hand, 
the structure could be undergoing alternative loadings , in which case either 
Q1 or Q2 would be acting but not together. Assume that the mode in which 
each component or element fails is known. The corresponding failure stress 
(of the ith element) is denoted by F^, which is taken to be normally dis- 
tributed with mean Fj_ and standard deviation cr-p.. It is required that the 

probability of failure of the structure not exceed a fixed value a. * 

Let the parameter Ai represent the size of a component to be deter- 
mined. The design parameters Ai constitute the variables for the optimiza- 
tion problem. For a given structure, the design stresses are functions of 
member sizes Ai and the elasticity E of the material. For the sake of 
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convenience, Ap and E are not taken to be random variables, but constants, 
in order to simplify the determination of the frequency distribution of the 
applied stress. The weight (W) of the structure is expressed as a function 
of the design parameters Ap . The optimization problem consists in determin- 
ing Ap so that W is a minimum and the probability of failure of the 
structure Pp is less than or equal to a. 

Let PFqi and ‘ p Fq 2 r ®P r esent the probabilities of failure of the 

structure under loadings Q1 and Q2, respectively. If the loadings Q1 and 
Q2 are such that they may or may not act simultaneously, then the structure 
can fail when either Q1 or Q2 or both Q1 and Q2 are acting on the struc- 
ture. Then, the desired probability of failure of the structure is the 
probability that it will fail under 

^FqI or Q2 or (Q1 and Q2) ^ an< ^ 


?F Q1 or Q2 or (Q1 and Q2) 


For any loading, the probability of failure of the structure (Pp) is 
expressed in terms of the probabilities of failure of its constituent 
elements (p^) as follows: 

pf = i - (i - pnu - p 2 ) . . . (1 - Pi) . . . (1 - p r ) (2) 


A structural element will fail when the stress (fp) in it, due to the 
applied loading, exceeds the allowable limit Fp. Let Pp denote the 
probability of failure of an element. Thus 


Pi = P(fi > Pi) = P(F i - fi < 0 ) ( 3 ) 


. or Q2 or Q1 and Q2 
m by 


acting together 


%1 + PF Q2 


- P Fqi P Fq2 


( 1 ) 


The quantity of interest is, then, the random variable (Fp - fp) . Let us 
consider the distribution of fp. It is a function of the loading and the 
stiffness of the structure. The latter quantity is assumed to be a constant, 
for the sake of simplicity. Hence, if the loading is normally distributed, 
fp will also be normally distributed. Also, Fp is normally distributed. 

The random variable (Fp - fp) is thus normally distributed. Now given the 
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values of the design parameters Aj_ , the random variables (F^ - f^) all become 
determinate; in other words , the means and standard deviations are known. It 
is then possible to determine the probabilities (PjO of the constituent 
elements and hence from equations (l) and (2) the probability of failure of 

the structure, that is. either P-p or P-m , as the case may be. 

*Q1 and Q2 ^Ql or Q2 

Thus, Pp is a function of the design parameters Aj_. 


P F “ 


P-ri or 

*Q1 and Q2 


' f Q 1 or Q2 or Q1 + Q2, 


J = 9(A 17 A 2 , 


A • , 


A r ) 


(4) 


The weight (w) of the structure can also be expressed in terms of the design 
parameters : 


W — W(A 1 ; A^>, • • • A j_ , ... AjO (5) 

The optimization problem consists in optimizing the function W (eq. (5)), 
subject to the constraint that Pp (eq. (4)) is less than or equal to the 
permissible limit a. Now that the optimizing function and the constraint 
have been expressed in terms of the design parameters, the stage is set for 
the optimizing process. To this end, the method of structural synthesis 
(ref. 2) can be employed. 


Optimization Method 

The structural optimization problem being considered is a mathematical 
programming problem for determining the optimum set of design parameters that 
minimizes the weight of the structure and restricts its probability of failure 
below a certain value. The method of structural synthesis is suggested for 
solving this problem. 

The synthesis method consists in repeated and systematic redesign of a 
structure by means of its governing principles without violating the con- 
straints imposed on the structure. A trial set of parameters is initially 
selected. For this set, the probability of failure of the structure under 
either Q1 or Q2 is calculated and checked to see if it exceeds the allowable 
limit a. If not, the assumed values of the design parameters are altered so 
as to diminish the weight. This modification is made according to the govern- 
ing principles of the method. With every modification, the probability of 
failure is calculated and checked as to whether it exceeds the allowable 
limit. The optimum set of parameters is said to have been reached when it is 
no longer possible to cause modifications in them, so as to diminish the 
weight of the structure without violating the probability of failure criterion 
P-p ^ cx. 

The method is qualitatively illustrated with a simple problem which 
requires only two design parameters to define a given design. Such a 
situation can be easily graphed so that the problem constraints and the 
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solution technique can be visualized* The design problem is to determine the 
optimum values of areas A x and A 2 of the two-bar truss (fig. l) so as to 
minimize (without the probability of failure exceeding a) the weight of the 
truss. Figure 2 is a plot in what may be called the design parameter space. 
The design parameters A x and A 2 are the abscissa and ordinate , so that any 
point in the design space represents a design. The weight contours 

W(Ai,A 2 ) = C m 

are plotted in the design space for different values (C m ) of the structural 
weight. The constraint curve 


Pp(Ai,A 2 ) - a 

is also shown plotted in the design space. Notice that the feasible region 
consisting of the set of all points of designs is bounded by the inequality 

P p (Ai,A 2 ) < a 

Clearly, the design satisfying the constraint which has the minimum weight 
occurs at M, the point of tangency of the minimum possible valued weight 
contour to the constraint curve. When there are more than two design param- 
eters, such a graphical solution is not possible. In such a case, the 
synthesis method provides means of systematic convergence to optimum design 
from a trial design. The two-dimensional case makes it possible to illustrate 
graphically the principles of this technique by means of figure 2. 

In terms of the design space, a design point could be 

(1) A free point (e.g., points ^ , n = 1,2, • - .), if it lies' 
in the feasible region 

(2) A bound point (e.g., point B) , if it lies on the constraint 
curve 

( 3 ) A violation point V, if it violates the constraint. 

In essence, the method consists of a combination of three directives in terms 
of effecting changes in the design parameters at a current point, depending on 
whether it is a free point, a bound point, or a violation point (see ref. 2 
for a more elaborate presentation of the method). They are as follows: 

(a) If the current design point is free (F n ) , changes are so made 
that the rate of change of weight is a maximum. This corresponds 
to a travel in the design space in the direction of the gradient 
of the weight function. 

(b) If the current design point is a bound point (B) , travel to an 
alternate point on the same weight contour and adopt directive 
(a). 
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(c) If the current design point is a violation (v), retreat to the 
previous design point and travel to an alternate point on the 
same weight contour from where it is possible to adopt directive 
(a) . 

Figure 2 shows a probable convergence path from a trial design point Fi to 
the optimum design point M. 

Optimization o f a truss for a given probab ility of failure constraint. - 
Consider a statically determinate truss AOB (fig* l) • The given data consist 
of (l) two loadings Q1 and Q2 which may or may not occur simultaneously. 

Both Q1 and Q2 are random variables , normally distributed with means Q1 
and Q2 and standard deviations Oqp and C7Q2.7 respectively; (2) failure 
stress Fj which is a random variable normally distributed with mean 7 and 
standard deviation tfpj and (3) the probability a that the structure will 
fail under either Q1 or Q2 or (Q1 and Q2), namely [PFqp Qr Q2 or ( Q1 and Q2 )]. 

The design parameters to be determined are areas Ai and A 2 of the bars 
OA and 0B> respectively. Areas A^. and A 2 are unknown constants. It is 
required to minimize the weight (w) of the truss subject to the constraint 
that the probability that it fails under Q1 or Q2 does not exceed a. 

If p is the density of the material and Z the length of each bar, 
then the weight (w) of the truss is given by 

¥ = p2(Ai + A 2 ) (6) 

Let Rxi and R 21 denote the forces in the bars OA and OB under the load 
Ql, and R 12 and R 22 , the forces in OA and OB under the loading Q2. Let 
P u a-nrl P 21 denote the probability of failure of bars OA and OB, respec- 
tively, under the loading Ql; and P 12 and P 22 denote the probabilities of 
failure of bars OA and OB, respectively, under the loading Q2. Row, 


P ij 



= P 





( 7 ) 


where i = 1, 2 for the bars OA and OB, respectively, and j = 1, 2 for the 
two load systems. Under the action of the load Q2, the bar OB is in com- 
pression and is, hence, liable to buckle. However, it is assumed that the bar 
OB has a low slenderness ratio and will consequently fail by yielding of the 
material at the stress F before the critical buckling stress is reached. 


Let g^j denote the random variable F - ( j/a ^ then 
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It is necessary to examine the random variable g. . to be able to determine 
P ,j_j. In general 


s u 



= F 


Ql 


Ci.(0) A i 


(9) 


where Cij(0) is a function of 0, Now since 9 is a constant, C-j_ j is a 
constant; also Aj_ is, in this distribution equation, a constant, and 
F and Qj are normally distributed. Hence, gj_j is a normally distributed 
random variable. The mean and standard deviation of the distribution g^j 
are given as mean 


3 ij 


= F 


c ij A i 


( 10 ) 


standard deviation 


a 


s u 



V v /g 

CifAi 2 / 


( 11 ) 


The constant C-j_ j is represented as follows: 


Cn = 2 cos 9 , 
C 2 i = 2 cos 9 , 
C 12 = 2 sin 9 , 
C 22 = 2 sin 9 , 

Now, the probability of failure 


for 

bar 

0A 

under loading 

Q1 

for 

bar 

OB 

under loading 

Q1 

for 

bar 

0A 

under loading 

Q2 

for 

bar 

OB 

under loading 

Q2 





<3g; 


1J 


The probability of failure of the structure is 


( 12 ) 


■^Ql or Q2 or (Q1 and Q2) 


= P Fq1 + PFq2 - Pf q1 PFq2 


(13) 


where 



(14) 

(15) 



(1 - P 13 .)(l - P 21 ) 
(1 - P 12 )(l - P 22 ) 


In the characteristics of the distribution 


g. . (eqs. (10) and (ll)), the 

X J 


quantities F, cr-p, Qj , Oq 4 , and are available from the given data. 


The 


A i* 


Thus, the characteristics of g^j are unknowns 


Hence , P* ^ 

x J 


only unknown quantity is 

in A that is, they are unknown functions in the variable 

(eq. ( 12) ) is also a function of the variable Aq. From equation ( 13 ), the 
probability of failure of the structure is an unknown function in the vari- 
ables Ax and A 2 (i.e., P-™ = Pp(A x , A 2 )). The condition 


Pp ~ BpUi, A 2 ) < a (l 6 ) 

is thus a constraint in the variables A x and A 2 . For the structural optimi- 
zation problem, the next function is the weight ¥ of the truss (given by 
eq. ( 6 )), the constraint is represented by equation (l 6 ), and the design 
parameters are A x and A 2 . The synthesis method can then be employed to 
obtain the optimum set of design parameters A x and A 2 so that W is a 
minimum and Pp does not exceed a. 


NUMERICAL EXAMPLES 


The above sample problem is solved for the following two sets of data: 

F = 15,000 lb/in. 2 , a F = 500 
Q 1 = 30,000 lb, Q 2 = 30,000 lb 
a Ql = a Q 2 ~ 1,000 
0 = 45 ° 

(1) P F = 0.01 (2) P F = 0.2 

The graphical solutions of the problem are presented in figure 3* The 
optimum design parameters for the two cases are: (l) Ai = 1.615 in. 2 , 

A 2 = 1.615 in . 2 with the weight parameter W/pZ = 3»23 in. 2 , and (2) 

Ai = 1.524 in. 2 , A 2 = 1*524 in . 2 with the weight parameter w/pZ = 3*048 in. 2 . 

On the other hand, if the idealistic version of failure (P F = 0) were 
adopted (i.e., if failure were prevented with certainty in designing the 
truss) , then extreme values of loading and allowable stress would have to be 
considered. Let Q F be the upper extreme value of loadings Q,1 and Q2 such 
that the probability of the load exceeding Qg is 0 . 987 X 10 " 9 ; then for the 
above data Qg = 36,000 lb. Let Fg represent the lower extreme value of the 
stress F such that the probability of the stress being less than Fg is 
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0. 987 x 10 " 9 ; then for the above data F E = 12,000 psi. For this case, then, 

the design parameters are Ai = 2.l4 in. 2 , A 2 = 2.l4 in. 2 , and the -weight 
parameter W/pZ = 4.28 in. 2 . 

CONCLUDING REMARKS 


The concept emphasized in this paper is the employment of probability of 
failure constraint in the mathematical programming problem to determine the 
set of design parameters which minimizes the weight of the structure. It is 
shown how the probability of failure constraint of a structure can be formu- 
lated in terms of the design parameters; then, knowing the optimizing weight 
function in terms of the design parameters defines the mathematical model of 
the structural optimization problem completely. 

It is seen that for the same loading, a truss designed for a probability 
of failure of 0.2 is lighter than that designed for a probability of failure 
of 0.01. However, both these structures are much lighter than that designed 
to prevent failure absolutely, that is, designed for Pp = 0. The probabil- 
ities of failure of 0.2 and 0.01 are realistic compared to the idealistic 
probability of failure of zero. If these probabilities are arrived at by a 
judicious compromise between the weight saving desired and the amount of risk 
that could be afforded, then it can be seen that considerable weight can be 
saved by designing the structure for realistic probabilities of failure. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 12, 1966 
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Area of member OA 



Figure 3* - Design parameter space showing the graphical solutions of the 
truss -optimization numerical examples. 
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